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ABSTRACT 

The  feasibility  of  utilizing  a  gaseous  core  nuclear 
reactor  to  provide  high  enthalpy,  high  pressure  gas  flow 
for  simulating  atmospheric  re-entry  conditions  was  inves- 
tigated.  The  test  facility  uses  a  mixture  of  nitrogen  and 
uranium  in  a  closed  cycle  with  no  attempt  to  contain  the 
uranium  fuel  within  the  core.   The  primary  purpose  of  the 
facility  is  to  provide  high  enthalpy,  high  shear  flows 
for  testing  re-entry  materials  and  shapes. 

Investigated  in  this  study  were  the  effects  of  the 
nitrogen-uranium  mixture  on  reactor  criticality,  nuclear 
contamination  of  the  test  model,  protection  of  the  reactor 
core  and  nozzle  structure  from  imposed  heat  loads  and  operating 
limitations  of  the  test  facility. 
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I     temperature,  Kelvin 
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density  proportionality  constant  ,  0.715 

1  3 
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t"  ratio  of  specific  heats 
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INTRODUCTION 

Many  problems  areas  have  been  uncovered  and  investi- 
gated in  our  attempt  to  master  space .   One  area  is  the 
re-entry  of  a  ballistic,  orbital  or  superorbital  body 
into  a  planetery  atmosphere  and  the  delivery  of  an 
payload  to  the  planet's  surface. 

The  kinetic  and  potential  energy  of  the  body  must 
be  dissipated  to  insure  successful  payload  delivery.   Atmos- 
pheric braking,  with  the  transfer  of  energy  from  the  high 

velocity  re-entry  body  to  the  surrounding  air,  is  often 

1 
the  optimum  entry  system. 

The  body  must  be  protected  from  the  high  heat  loads 
due  to  friction  and  the  rapid  compression  heating  of  the  air 
immediately  surrounding  the  body.   Ablative  materials  are 
used  to  both  protect  the  body,  by  acting  as  an  insulator, 

and  to  carry  away  the  heat  that  is  absorbed  by  the  ejection 

2 

of  mass  of  the  ablative  material. 

The  ablative  proqess,  however,  is  very  complex  and 

attempts  to  analytically  predict  the  behavior  of  new 

3 

materials  and  shapes  have  met  with  little  success.   Various 

attempts  at  simulating  re-entry  conditions  have  been  made, 
but  no  existing  facility  can  completely  duplicate  the 


re-entry  environment,  although  some  aspects  of  the  problem 

4,  5 

can  be  simulated. 


Existing  test  facilities  cannot  duplicate  the  high 
enthalpy,  high  shear  regime  for  a  sufficiently  long  test 
time  to  be  able  to  investigate  this  critical  ablation  region. 

This  study  will  investigate  the  feasibility  of  a  high 
enthalpy  test  facility  which  uses  a  fissionable  gas  reactor 
as  a  power  source.   The  gas  core  reactor  drives  a  closed 
cycle  test  system  and  provides  the  necessary  energy  for 
high  enthalpy,  high  pressure  gas  flows  of  long  duration. 

In  establishing  the  feasibility  of  this  system,  the 
re-entry  problem  is  first  investigated  to  establish  the 
power  requirements  of  the  system. 

Stagnation  pressure  and  enthalpy  are  determined  by 
the  method  outlined  in  Reference   1,   from  these  the 
power  required  to  simulate  an  atmospheric  entry  in  a  ground 
test  facility  is  determined. 

Then  the  reactor  itself  is  studied.   The  gaseous  core 

reactor  utilized  is  similar  in  operating  principle  to 

the  gaseous  core  nuclear  rocket  engine  proposed  as  the 

6,7,8 
next  generation  space  booster.       However,  as  the 

proposed  system  incorporates  a  closed  operating  cycle, 

containment  of  the  nuclear  fuel  does  not  present  an  obstacle 

as  it  does  in  the  gaseous  core  nuclear  rocket.   Critical 

size,  fuel  concentration  and  fuel-nitrogen  ratio,  where 

nitrogen  is  the  working  fluid,  are  determined  using  the 


9 

diffusion  approximation  and  one  group  neutron  theory. 

The  contamination  problems'  associated  with  the  utilization 
of  a  nuclear  medium  are  also  investigated. 


The  combined  heat  loads  to  the  reactor  and  nozzle 

10 

are  examined.   The  heat  loads  in  the  nozzle  are  determined 

and  homogenous  transpiration  cooling    to  protect  the  nozzle 
is  examined.   The  use  of  a  seeded  carbon  particle  layer 

to  protect  the  moderator  from  radiation  heat  loads  is  also 

.   ,  12,  13 
examined. 


Finally  the  power  requirements  of  auxiliary  equipment 
to  maintain  steady  state  closed  cycle  operation  are  examined 


CHAPTER  I 


RE-ENTRY  CONDITIONS 


1.1   Pressure  and  Enthalpy 

To  test  the  materials  used  in  a  re-entry  vehicle  it 
is  necessary  to  know  the  velocity,  the  pressure  and  the 
heating  involved  in  the  re-entry. 

As  a  typical  example  consider  a  vehicle  returning 
to  earth  with  a  superorbital  velocity.  An  estimate   of 
the  velocity  can  be  obtained  by  assuming  it  equal  to  the 
velocity  needed  to  escape  from  a  nonrotating  earth,  with 
no  drag  forces.   The  energy  per  unit  mass  to  escape  from 
the  earth,  acting  as  a  central  force  field  is 

■^\4a    =     -£  d.i) 

and  the  velocity  to  escape  is 


*  --  A 


&r 


(1.2) 


when  k   is  the  earth  gravitational  constant,  A*\o   ""Vsec* 
and   R.  is  the  earth's  radius  ,  <©.  4.  x  lo&  ro.  ,    the  escape 


velocity  from  the  earth's  surface  is  1.12  x  10   m/sec 


This  is  the  minimum  for  earth  escape  and  is  the  mini- 
mum at  which  a  superorbital  body  will  approach  the  sensible 
atmosphere,  commencing  re-entry. 

Expressions  for  stagnation  pressure  and  enthalpy,  at 
the  point  of  maximum  pressure  and  also  at-  the  point  of 
maximum  heating  rate,  may  be  found  by  using  the  velocity 
altitude  history  given  in  Reference  l  for  an  exponential 
atmosphere, 


(1.3) 


v  «  vc  e 


t  -cP  a  Po  «  g   :  i 


where  Ve  is  the  initial  entry  velocity,  Co  is  the  coeffi- 
cient of  drag  of  the  body,   m   the  mass  of  the  body,   ©e 
the  entry  angle  and  %      and  fi    are  constants  defined  by 

-**  (1.4) 


f  =  f - 


<k  e 


The  stagnation  pressure  on  the  body  is  approximated 
by  Newtonian  theory  with  the  body  angle  considered  small 
and  the  ambient  pressure  neglected.   The  stagnation  pres- 
sure at  any  point  in  the  entry  is 


-A 


-Aiwa   f-C»K  P°«e"^  N 


For  maximum  stagnation  pressure,  set    -r -  O    and 


(1-5) 


obtain 


^ 


£  - Qe-^)eC*e"*\  C.e  V^ C»e*Mft)  -O    (i.6) 

where 

G  h  ^(XN/*1 


Cv   and    Ci.    are  fixed  by  the  atmospheric  model 
ballistic  parameter,  m|>,  ge' an<^  initial  entry  velocity; 
and  remain  constant  during  re-entry. 

The  altitude  for  maximum  stagnation  pressure  during 
re-entry  is  given  by 


i-ii*-;^    ¥■         a.,) 


Since  the  velocity  for  maximum  stagnation  pressure  is 
given  by 


V  -Ve  e 


*  (1.8) 


and  associated  density  is 


the  maximum  stagnation  pressure  is 


(1.9) 


o 


***        Co  A, 


Sxo  $e    AJye  (1.10) 


7 
and  the  stagnation  enthalpy  is 


Mo  =  la  *  -5-  (l.lla) 


Z 

or  since  the  static  enthaply  is  very  much  less  than  dynamic 
enthalpy, 

U  =  ¥±-  (l.llb) 

Similarly,  expressions  are  obtained  for  stagnation  pressure 
and  enthalpy  at  the  point  of  maximum  rate  of  convection 
heating  of  the  stagnation  point,  where  stagnation  point 
heating  rate  is  given  by 


where  C^  is  the  stagnation  point  coefficient  of  heat  trans- 
fer, a  function  of  body  configuration,   <T"   is  the  radius  of 
curvature  at  the  stagnation  point,  and   P   and  V  are  functions 
of  altitude.   To  obtain  the  point  of  maximum  stagnation  point 
heating  rate,  the  previously  derived  expressions  for   P 
and   V    are  substituted  in  the  above. 


d    /dU.\  „ 

For  the   maximum   set    —-  [  ^  ')sOpbt 


cl 


<JLV  L»»      V  d0  k         3ecr  V« 


c    (l.ia) 


which  occurs  at  an  altitude  of 


(1.14) 
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and   velocity   of 


-7<« 

v  -  vc    e  (i.i5) 


The  associated  stagnation  pressure  is 


O  -  ms^Ge  Ve*  A  (1.16) 


and  stagnation  enthalpy  is  given  by 


Ho  ■  ^V  (1*17) 


Figure  i  shows  the  extremes  of  stagnation  pressure 

verses  stagnation  enthalpy  which  are  encountered  during 

an  atmospheric  re-entry.   The  ballistic  parameter,  JS  Sift 6c    / 

CoAk 

is  a  free  parameter  in  determining  the  envelope.   Two 

4 

envelopes  are  shown,  one  for  Ve*  1.12*  10  to/sec   corresponding 
to  velocity  entering  the  sensible  atmosphere  equal  to 
escape  velocity.   The  other  is  for  Ve  -  O.S  *  lO  y^/sec, 
corresponding  roughly  to  lOo  r\.r\.      orbital  velocity.   It 
can  be  readily  seen  in  Figure  1    that: 

i)  stagnation  enthalpy  at  maximum  heating  rate  and 
maximum  stagnation  pressure  is  independant  of 
the  ballistic  parameter,  but  increases  with  in- 
creasing initial  entry  velocity. 

ii) stagnation  pressure  increases  directly  as  entry 
velocity  and  directly  with  ballistic  parameter. 

1.2  Power  Requirements  for  Simulation 

The  stagnation  pressure  on  the  test  model  is  fixed 


by  re-entry  conditions  but  the  chamber  pressure  is  limited 
by  structural  considerations.   The  maximum  chamber  pressure 
is  set  at  1000  atmospheres  for  the  analysis.   It  is  desir- 
able to  test  at  the  lowest  possible  Mach  numbers  consistent 
with  representing  hypersonic  pressure  and  heating  rate  dis- 
tributions.  For  this  reason,  and  to  minimize  nozzle  losses 
and  cooling  requirements,  a  test  Mach  number,  Ht  ,  of  3.0 
is  used  in  this  analysis. 

For  isentropic  frozen  flow  through  the  nozzle,  the 
power  required  per  unit  throat  area  is 


p         p.h.       / —       -  uvl — ■ 


If  Equations  1.16  and  1.17  are  used  for    "v^     and 
^0   respectively,  the  power  required  per  unit  throat  area 
to  simulate  the  conditions  at  the  point  of  maximum  heating 
rate  is 


A»     /2KT?   c 


& 


These  power  requirements  are  plotted  in  Figure  2  against 
initial  entry  velocity  with  the  ballistic  parameter  as  a 
free  parameter.   Reference  14  was  used  to  obtain  the  charac- 
teristics of  nitrogen, which  was  taken  as  the  test  gas. 

Limiting  the  chamber  pressure  to  1000  atmospheres 
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established  one  limit  on  the  re-entry  profiles  that  can 
be  simulated.   Another  is- set  by  limiting  the  chamber 
temperature  to  10,000  K,  this  upper  limit  is  set  after 
a  consideration  of  the  heat  transfer  loads  experienced 
in  the  various  test  facility  components. 

1.3  Test  Time 

The  main  purpose  of  this  facility  is  to  examine  the 
behavior  of  re-entry  materials  and  shapes  under  the  high 
heat  and  shear  loads  of  an  atmospheric  re-entry.   Such 
materials  are  ablative  in  nature,  that  is,  they  absorb 
heat  by  increasing  in  temperature  and  changing  in  chemical 
or  physical  state.   The  heat  is  carried  away  from  the 
surface  by  a  loss  of  mass.   The  departing  mass  also  blocks 

part  of  the  convective  heat  transfer  to  the  remaining  mate- 

2 

rial.   The  tunnel  running  times  are  set  by  the  time  required 

for  the  rate  of  ablation  to  reach  equilibrium.   This  time 

15 
is  a  function  of  the  ablative  material    and  for  typical 

2 
ablative  materials   ,  such  as,  graphite  and  glass  reinforced 

phenolics,  the  ablative  equilibium  times  is  extended 

to  200  seconds . 

1.4  Existing  Facilities 

Shock  tubes  are  used  to  generate  higher  temperatures 

and  pressures  to  simulate  re-entry  conditions.   The  Cornell 

16 
Aero  Laboratory  six  foot  hypersonic  tunnel  can  generate 

o 

stagnation  temperatures  and  pressures  up  to  4200  K  and 

2000  atmospheres  respectively.   The  test  time  for  these 
pressures  however  is  only  4  micro-seconds.   This  tunnel  can 
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be  upgraded  to  10,000  K  and  2000  atmospheres  but  still  has 
a  run  time  far  too  short  to  produce  the  time  history  of 
high  temperatures  and  pressures  at  which  the  ablative 
materials  fail. 


The  test  time  for  a  shock  tube  may  be  increased  by 
using  several  tubes  fired  in  a  Gat ling  gun  fashion.   This 
however  only  increases  run  time  to  seconds,  developing 
a  jet  power  of  25,000  HP  ,    about  four  orders  of  magnitude 
too  small  for  re-entry  requirements. 


High  enthalpy  and  pressure  flow  is  also  possible 
utilizing  a  hot  shot  tunnel  powered  by  an  electrical 
discharge  arc.   This  type  of  tunnel  requires  large  amounts 
of  stored  electrical  energy  and  is  presently  limited  to 
a  few  hundredths  of  a  second  operation. 

Although  these  facilities  can  approach  the  desired 
levels  of  enthalpy  and  pressure,  they  cannot  maintain  these 
conditions  for  a  sufficient  lenght  of  time.   They  are  severely 
power  limited.   As  an  example,  the  power  required,  from 

Figure  2,  to  simulate  the  re-entry  of  a  vehicle  whose  ball- 

2 

istic  parameter  is  2000  grams/  cm     and  initial  entry 

/  2 

velocity  is  8000  meters/second,  is   30   megawatts  per  cm  of 

throat  area.   This  is  about  ten  percent  of  the  total  output 
of  the  Grand  Coulee  Dam  for  each  square  centimeter  of  nozzle 
throat  area. 

Plasma  arc  tunnels  provide  high  enthalpy  but  low 
pressure  flows  and  are  currently  being  used  to  test  re-entry 
materials.   These  facilities  are  also  power  limited  and  operate 
at  low  stagnation  pressures  and  use  argon  as  a  working  fluid 
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to  achieve  long  duration  tests  with  comparatively  low 
powe  r  input  s  . 


Because  of  these  high  power  requirements,  a  gaseous 
core  reactor  is  suggested  to  drive  the  test  facility-/  the 
gaseous  core  because  its  operating  temperature  is  limited 

only  by  heat  loads  to  the  containing  structure  and  can 

o 

operate  at  the  desired  temperature,  10,000  K. 
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CHAPTER  II 
TEST  FACILITY 

2.1  Gaseous  Core  Reactor  Requirement 

The  power  required  to  produce  long  duration  high 
pressure ,high  enthalpy  gas  flows  appears  to  he   technically 
achievable  only  with  nuclear  power.   Conventional  combus- 
tion is  limited  by  the  recombination  of  combustion  products , 
which  absorb  energy  and  result  in  an  equlibrium  tempera- 
ture far  below  the  required  stagnation  temperature.   In 
atomic  reactors,  solid  and  liquid  core  operating  tempera- 
tures have  upper  limits,   the  former  is  the  melting  point 
of  the  nuclear  core  and,  or,  structure  and  the  later  is 
boiling  point  of  the  liquidized  fuel,  and,  or, melting  point 
of  the  structure. 

A  gaseous  core  nuclear  reactor  is  capable  of  steady 
state  operation  at  very  high   temperatures  and  pressures. 
The  operating  temperature  of  a  nuclear  reactor  utilizing 
a  fissionable  gas  core  is  limited  only  by  the  heat  load 
that  the  hot  gas  imposes  on  the  container.  These  heat  loads 
can  be  controlled  to  some  degree  by  a  combination  of  con- 
vective  and  transpiration  cooling  and  by  injection  of  a 
thermal  radiation  absorber  to  block  radiation  heat  transfer 
to  the  wall.   The  gaseous  core  reactor  in  this  study  is 
capable  of  providing  the  energy  necessary  to  run  a  test 
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at  the  desired  stagnation  pressures  and  enthalpy,  that  is, 
7,100  megawatts. 

2.2  Gas  Cycle 

To  simulate  the  effects  of  atmospheric  heating,  nitrogen 

is  used  as  the  working  medium.   Uranium  235,  initially  a 

fine  power,  is  the  nuclear  fuel.   Biological  hazards 

and  economic  necessity  make  it  mandatory  that  the  uranium 

fuel  be  contained  within  the  core  of  a  gas  core  nuclear 

18 
rocket.    As  the  proposed  system  is  a  closed  system,  these 

considerations  do  not  apply.   This  is  not  to  say,  a  priori, 

that  fuel  containment  in  the  core  proper  is  not  required. 

Presence  of  uranium  in  the  nitrogen  flow  complicates  the 

analysis  and  interpretation  of  experimental  data.   If  the 

ratio  of  mass  flow  of  uranium  to  nitrogen  is  limited  to 

a  maximum  value  of  ten  percent,  or  about  one  percent  mole 

ratio,  it  should  permit  the  assumption  that  the  flow  in 

the  nozzle  and  test  section  behave  as  a  pure  nitrogen  flow. 

This  limit  is  imposed  as  the. maximum  allowable  uranium 

content  for  thermodynamically  simulating  the  effects  of 

atmospheric  entry. 

The  presence  of  fissionable  gas  in  the  closed  cycle 
causes  induced  radioactivity  of  the  system  components. 

State  of  the  art  shielding  and  radioactive  materials  handling 

9 

techniques  are  assumed  sufficient  to  provide  necessary 

biological  protection.   As  this  facility  would  be  very 
expensive  to  construct, it  is  necessary  that  the  induced 
radioactivity  of  the  model  be  within  biological  limits 
for  examination  of  the  model  within  a  short  time  after 
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the  test  and  associated  radioactive  exposure. 

The  radioactive  damage  or  effects  to  the  model  must 
be  slight  enough  so  that  these  associated  effects  do  not 
mask  or  override  the  aerodynamic  effects  which  are  being 
investigated . 

These  nuclear  effects  on  the  model  are  investigated 
to  determine  their  influence  on  the  design  of  the  over-all 
system,  and  to  determine  if  containment,  greater  than  that 
implied  by  the  thermodynamic  limit  of  ten  percent  uranium, 
is  required . 

Chamber  pressure  and  temperature  are  set  at  1000 

o 

atmospheres  and  10,000  K.   A  portion  of  the  heat  load  to 

the  moderator  is  absorbed  by  the  flow  of  nitrogen  entering 
the  core,  the  remainder  is  removed  by  an  external  water  flow 
enroute  to  another  system  component,  the  water  diffusor.   A 
seeded  layer  of  carbon  is  used  to  provide  a  radiation  barrier 
which  protects  the  moderator  from  radiation  heat  loads. 

The  thermal  energy  of  the  hot  gas  in  the  reactor  core 
is  converted  to  kinetic  energy  in  a  convergent-divergent 
nozzle.   The  nozzle  is  cooled  by  conventional  backside 
cooling  and  transpiration  cooling,  both  utilizing  a  secondary 
flow  of  high  pressure,  but  low  enthalpy  nitrogen. 

A  schematic  of  the  test  facility  is  given  in  Figure  3. 
A  temperature  -  enthalpy  diagram  of  the  nitrogen  cycle  is 
given  in  Figure  4. 
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After  the  test  section,  the  secondary  nitrogen  flow 
is  vented  into  the  main  flow,  reducing  the  temperature 


slightly.   The  stagnation  temperatures  are  still  in  the 

o 

vicinity  of  10,000  K  and  prevent  the  use  of  conventional 

turbo-expansion  devices  or  heat  exchangers  at  this  point 


A  water  spray  diffusor  is  utilized  to  slow  the  flow 
and  lower  the  stagnation  pressure  and  temperature  of  flow 
to  1000  K  and  100  atmospheres.   The  water  spray  is  sufficient 
to  cool  the  walls  of  the  diffusor. 


A  centrifugal  separator  is  used  to  separate  the 
condensed  uranium  particles  from  the  flow  of  superheated 
steam  and  nitrogen  leaving  the  water  diffusor. 

At  this  point  the  flow  is  basically  free  of  uranium 
and  at  a  temperature  and  pressure  at  which  conventional 
turbo-machinery  can  be  employed.   The  flow  has  2,500  megawatts 
of  available  power,  if  expanded  only  to  the  saturation  point 
of  the  water.   After  this  energy  is  removed  from  the  stream, 

the  water  is  condensed  and  separated.   The  nitrogen,  at 

o 
0.24  atmospheres  and  370  K,is  cooled  further  in  a  conven- 
tional heat  exchanger  and  pumped  back  to  1000  atmospheres 

using  two  multistage  centrifugal  flow  compressors  with 

o 

interpump  cooling  to  maintain  temperatures  below  1000  K. 

The  power  required  for  the  pump  work,  880  megawatts,  is 
provided  by  the   aforementioned  turbine. 

The  pressure  of  the  water  is  raised  to  120  atmospheres 
in  another  multistage  centrifugal  pump,  requiring  10  mega- 
watts of  the  above  turbine  work. 
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The  temperature  -  enthalpy  diagram  of  the  water  cycle 
is  given  in  Figure  5. 

2.3   Component  Interface  Problems 

To  be  successful  as  a  system,  the  components  of  the 
facility  must  be  mutually  compatible.   In  the  following 
sections,  the  operating  characteristics  of  the  reactor 
and  the  nozzle,  which  are  the  main  components  of  the 
system,  are  examined.   The  reactor  is  examined  to  determine 
the  effects  of  various  operating  pressures  and  fuel-nitrogen 
concentration  on  the  critical  size  and  fuel  concentration. 

Contamination  of  the  model  is  examined  to  determine 
if  it  imposes  a  constraint  on  the  fuel-nitrogen  mixture 
in  the  core . 

Heat  transfer  in  the  reactor  and  nozzle  are  also 
examined  to  determine  limits  on  operating  temperature, 
power  and  component  size. 
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CHAPTER  III 


GASEOUS  CORE  NEUTRONICS 


3.1  General 

Reference  6  shows  that  the  gaseous  core  nuclear 
reactor  is  theoretically  feasible.   It  shows  that  one 
can  obtain  a  reasonable  estimate  of  reactor  critical 
radius  and  critical  fuel  concentration  by  considering 
the  reactor  core  to  be  a  homogeneous  mixture  of  fuel 
and  moderator  or  a  reflected  sphere  of  fuel  alone.   The 
effects  of  containing  the  fuel  with  a  gaseous  annulus 
of  hydrogen  are  also  investigated  as  were  the  effects  of 
multiple  cavities.   Reference  6  quantitatively  shows  that 
a  gaseous  core  reactor  can  be  designed  with  either  simple 
or  multiple  cavity  geometry. 

In  this  study,  the  reactor  will  be  considered  to  be 
a  homogeneous  mixture  of  pure  uranium  23  5  and  Ni  with 
a  bare  spherical  geometry.   The  effects  of  reflecting 
the  above  bare  sphere  are  also  considered.   One  group 
neutron  theory  is  utilized,  that  is,  the  neutrons  are 
considered  to  be  in  thermal  equilibrium  with  the  reflector- 
moderator. 

No  attempt  is  made  to  separate  the  fuel  from  the 
working  fluid,  nitrogen.   Criticality  calculations 
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are  made  with  the  ratio  of  number  density  of  nitrogen 
atoms  to  number  density  of  fuel  atoms,  NIm  /Mu  t 
as  a  parameter.   This  gives  one  measure  of  contamination 
of  the  flow  in  the  model  test  section.   From  the  fluid 
dynamic  point  of  view,  the  flow  is  treated  as  behaving 
as  pure  nitrogen,  and  the  number  density  ratio  gives  a 
measure  of  how  valid  this  approximation  is.   The  model 
is  also  subjected  to  radioactive  contamination,  the 
number  density  ratio  gives  a  partial  measure  of  this  effect; 
partial  in  that  the  uranium   itself  in  the  flow  is  not 
the  only  source  of  model  contamination. 

3.2   Criticality  Condition  , 

As  a  first  approximation,  the  reactor  can  be  considered 
to  be  a  bare  homogenous  mixture  of   pure  uranium  23  5  and 
nitrogen.   Criticality  is  determined  by  considering  all 
of  the  neutrons  in  the  core  to  be  thermal,  that  is,  there 
is  no  loss  of  neutrons  in  the  process  of  slowing  down 
from  the  high  energy  state  of  production  in  the  fission 
process,  to  capture  at  a  lower  thermal  energy  state.   For 
steady  state  operation,  the  conservation  of  neutrons  requires 
that  the  rate  of  production  of  neutrons  in  a  control  volume 
be  equal  to  the  rate  of  loss  of  neutrons  due  to   absorption 
and  leakage  from  the  volume. 

The  rate  of  production  of  neutrons  per  unit  volume 
is  given  by 


Production    =  V*  T  Zx+ 


3.1 
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where    N^    is  the  average  number  of  neutrons 

produced  per  thermal  fission 
&        is  the  neutron  flux 
and     ^f\t   is  the  macroscopic  fission  cross  section 

of  the  core  material. 

The  neutron  leakage  is  given  by  the  diffusion  approxi- 
mation as 

Leakage    -  -  D  ^  $  (3.2) 


where    V    is  the  La  Placian  operator 
and      D    is  the  diffusion  coefficient, 

1/3  of  the  transport  mean  free  path 

in  the  reactor 

The  absorption  term  is  similar  to  the  fission  term, 


Ab 


sorption  -     ^J  (  £  a,  +£*V  (3.3) 


here     2jx      **  the  macroscopic  absorption  cross  section 
in  the  core . 

Equations  3 .1, 3 .2  and  3.3  are  combined  and 
the  conservation  of  neutrons  is  expressed  in  differential 
form  by 

DvV  *  (v>£f-£oJ^   -  O  (3.4) 

or  if  the  material  buckling,  8M  , is  defined  as 

a,2  =  Jitt-ta.)/    D  (3.5) 
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the  conservation  of  neutrons  becomes 


v*<^  +  B2,  <j6  *  o 


(3.6) 


For  a  spherical  reactor,  with  the  boundary  conditions 
that  the  flux  goes  to  zero  at  the  radius  of  the  core  and 
the  gradient  of  the  flux  is  zero  at  the  center  of  the 
sphere.   Equation  3.6   can  be  solved  for  the  geometric 
buckling  of  the  reactor, 


a.1  -  «a 


(3.7) 


1  VR> 

The  criticality  condition  is  the  matching  of  the 
geometric  buckling  which  satisfies  Equation   3.6    with 
the  material  buckling  of  Equation   3.5.        The  diffusion 
length  is  defined  as 


(3.8) 


and  the  infinite  multiplication  factor  as 


With  this  notation,  the  criticality  requirement  is 


w*v 


(3.10) 


The  criticality  requirement  for  the  system  is  that 
the  geometric  buckling  be  equal  to  the  material  buckling. 
This  means  that  the  rate  of  neutron  creation  due  to  fissions 
is  equal  to  the  rate  of  destruction  of  neutrons  due  to 
absorption,  that  is,  fission  and   nonfission  capture,  and 
leakage  from  the  critical  region. 


For  a  nitrogen  and  pure  uranium  235  system,  the  number 


of  neutrons  released  per  thermal  fission  of  uranium  235  is 
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MS 


The  interaction  cross  sections  per  atom  for  uranium, 


sub-subscript   u.  /  and  nitrogen,  sub-subscript   •* 


are 


18 


CkyK  =  580  barns 
c-  =  680  barns 
Qr   =1.88  barns 


Gv*.  =  10  barns 
C4m  =  10  barns 


where  the  subscripts 


■f.   denotes  fission, 
ol.   denotes  absorption, 
g   denotes  scatter. 


With  these  cross  sections,  the  criticality  equation 
can  be  used  to  determine  the  number  density  ratio  of 
nitrogen  atoms  to  uranium  atoms  required  for  criticality 
as  a  function  of  the  core  radius  and  the  number  density 
of  uranium  235  atoms, 


fe  '  199  -A 


4 
*10  - 


ISAS  x  ID 
R1  Ma* 


47 


(3.11) 


Using  this  relation  the  critical  radius  is  found  as 
a  function  of  uranium  235  concentration  with  Mm  /r4u 
as  a  parameter.   The  results  are  given  in  Figure  6  with 
the  mass  flow  ratio  of  nitrogen  to  uranium  shown  in  lieu 
of  the  number  density  ratio. 


Examination  of  the  criticality  equation  above  shows 
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that  as  Mn /Mo  decreases  below  199,  R  si    must  increase. 
Therefore,  fsj*  /\l,j    =  199  gives  the  minimum  reactor  size 
for  a  given  concentration,  and  the  associated  mass  flow 
ratio  is  15.4.   This  is  due  to  the  scattering  contribu- 
tions of  the  nitrogen  atoms  in  the  core  reducing  the 
transport  mean  free  path  of  the  neutrons.   This  is  most 
readily  seen  if  the  criticality  condition  is  expressed  by 

0*  _    IT2  3(1-  <io<e} 

*   ~    Mo1  (  v>  <r,     -  ^  -  ^  <w  \(  <r*o  *  *£m«  VZ  )        (3-12) 

where  the  macroscopic  cross  sections  are  replaced  by 
the  number  density  times  the  microscopic  cross  section. 
The  denominator  has  a  maximum  value  when  KM  /  slu      is 
199.   For  higher  or  lower  number  density  ratios,  the 
critical  radius  increases. 

As  used  here  R  is  not  the  actual  radius  of  the 
homogenous  sphere,  but  rather  that  radius  plus  the  extrapolation 
distance.   The  extrapolation  distance   is  the  distance 
at  which  the  neutron  flux,  which  has  a  finite  value  at 
the  physical  boundary  of  the  sphere,  goes  to  zero.   Reference 
9  gives  the  extrapolation  distance  as 

^r  =  2.\"5D  (3.i3) 

The  physical  radius  of  the  bare  sphere  is 

r  -  R-  2.13  0  (3*14) 

This  critical  size  and  uranium  concentration  so 
obtained  is  only  an  approximation  and  is  valid  only  in 
determining  the  effect  of  nitrogen  poisoning  of  the  reactor 
on  critical  size  and  critical  uranium  concentration.   The 
actual  reactor  will  be  operating  at  elevated  temperatures 
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which  will  effect  the  interaction  cross  section, which 
generally  vary  inversely  with  neutron  energy.   A  higher 
equilibruim  temperature  means  higher  neutron  energies 
and,  in  general,  lower  interaction  cross  sections.   Also 
the  actual  reactor  will  not  be  bare  but  must  be  enclosed 
in  some  way. 

The  reactor  core   is  in  fact  contained  in  a  pressure 
vessel, for  structural  reasons  the  walls  of  which  must  be 
maintained  at  a  low  temperature  relative  to  the  reactor 
core.   Lining  the  pressure  vessel  with  beryllium,  a 
moderator,  which  must  also  be  cooled  for  structural  reasons, 
provides  a  means  of  decreasing  leakage  from  the  core  and 
a  means  of  providing  low  energy  neutrons  to  the  core. 

Nitrogen,  being  heavier  than  beryllium,  does  not 
absorb  as  much  of  the  neutrons  energy  in  a  collision, 
therefore  more  nitrogen  collisions  are  necessary  to  thermalize 
a  neutron.   Also  nitrogen,  even  at  very  high  pressures,   is 
still  very  much  less  dense  than  the  solid  beryllium  moderator 
and  at  the  higher  core  temperatures,  the  scattering  cross 
section  for  nitrogen  is  less  than  that  for  the  cold  beryllium 
moderator.   It  is  reasonable  to  assume  then,  that  the 
neutrons  are  intharmal  equilibruim  with  the  cool  reflector 
and  not  with  the  hotter  nitrogen-uranium  mixture. 

The  effects  of  the  reflector-moderator  on  core  radius 
can  be  found  from  a  one  velocity,   two  region  consideration. 
As  before  in  the  core,  the  conservation  of  neutrons  is 
given  by  Equation  3.6. 
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In  the  reflector,  there  is  no  fissionable  material 
and  the  conservation  of  neutrons  gives 


Dr  sr*<f>r   -  fLr  $r     =  o 


(3.15) 


where  subscripts  r  and  c  stand  for  reflector  and  core 
respectively . 

With  the  inverse  of  the  square  of  the  reflector 
diffusion  length,  l_ir/  defined  as 


Kr      =  ***-*/ Or  (3-16) 


the  neutron  conservation  equations  in  the  core  and 
reflector  is  obtained  by  dividing  Equation  3.15  by  Dr-  , 

<?*&     *  Be  #c   -o  (3#17) 

and 

V*^r    ~  k.r    <fir      sO  (3.18) 


This  set  of  differential  equations  can  be  solved 
for  a  criticality  condition  with  the  boundary  conditions 

O  <  ^cCr-o)    <co  (3-.19a) 

#clr*Rr)    -    ^r(r-RO  (3.19b) 


***Uu  -    *r*L-* 


(3.19c) 


and 


^r(Rr*dc  +  S 1   "  O  (3.19d) 
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where   d*.     is  the  reflector  thickness 

and     R,-     is  the  critical  radius  of  the  reflected  core. 

The  solution  of  Equation   3.17  and   3.18   with 
boundary  conditions  of  3.19   gives  an  expression  for 
the  criticality  of  a  reflected  spherical  reactor  in  the 
transcendental  form, 

a(krR^U>    k.Uc    +*)"*    l)  vDc(l-RrBc    CoV&cRc)  O.20) 


The  reflected  core  radius  can  be  expressed  in  terms 
of  the  bare  core  radius,    r      ,        and  the  reflector  saving,  S, 
as 

R^  *  r"S  (3.21) 

This  value  is  substituted  into  Equation  3.20      and 
a  transcendental  expression  for  reflector  saving  is 
obtained, 

Col  US  =  Or  U   co^Vx  (<Ar  *3^  -  0-jg)  (3.22) 

Equation  3.22  is  solved  graphically  for  the  reflector 
saving  at  each  known  r,  then  the  reflected  core  radius,  R,.,  is 
obtained  by  Equation  3.21 

The  resulting  critical  radii   are  plotted  in  Figure  .7 
against  critical  uranium  concentrations  for  a  100cm 
beryllium  reflector.   Shown  is  the  reflected  core  radius, 
with  the  mass  flow  ratio  of  nitrogen  to  uranium  as  a  free 
parameter.   Also  shown  are  lines  of  constant  pressure  for 
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an  operating  temperature  of  10,000  K. 


Heat  transfer  in  the  reactor  core,  coupled  with 
practical  limits  on  the  capacities  of  the  system's  pumps 
and  turbines,  lead  to  the  choice  of  a  minimum  core  radius. 

Structural  consideration  dictates  a  low  operating 
pressure,  but  this  is  not  compatible  with  the  desire 
to  simulate  the  high  stagnation  pressures  of  a  re-entry. 
Also,  the  problems  associated  with  heat  transfer  in  the 
core  require  high  mass  flows  and  low  temperatures.   For 
these  reasons  the  operating  pressure  is  set  at  the  maximum 
possible  with  projected  structural  technology,  that  is, 
1000  atmospheres. 

While  it  is  also  desired  that  the  mass  flow  rate  of 
nitrogen  to  uranium  be  maximized  so  that  the  test  facility 
simulate  atmospheric  conditions,  the  requirements  for  a 
minimum  core  size  dominate  in  this  case.   It  is  also  noted 
that  the  critical  core  radius  increases  rapidly  with  a 
slight  increase  in  mass  flow  ratio  from  the  minimum  of  15.4 

For  these  reasons  the  core  radius  is  set  at  the  minimum 

possible  with  a  pressure  limit  of  1000  atmospheres.   This 

radius  is  170  centimeters,  occuring  at  the  mass  flow  ratio 
of  15.4 
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CHAPTER  IV 


MODEL  NUCLEAR  CONTAMINATION 


4.1   Radiation  Damage 

In  the  gaseous  core  nuclear  reactor  used  as  a  power 
source  for  this  system,  there  is  no  provision  made  in 
the  core  itself  for  the  separation  of  fuel  and  fission 
fragments  from  the  working  fluid.   This  makes  the  design 
of  the  reactor  much  more  simple  then  if  some  hydrodynamic 
means  of  separation  of  fuel  from  the  working  fluid  is 
utilized.   As  the  system  is  closed  there  is  no  problem 
of  venting  radioactive  material  to  the  atmosphere  Relieving 
one  of  the  major  constraints  that  faces  the  use  of  a  gas 
core  nuclear  rocket. 

The  presence  of  radioactive  material  in  the  fluid 
stream  causes  contamination  of  the  model  being  tested. 
This  contamination  is  examined  to  determine  if  the  structural 
damage  to  the  model  will  signicantly  influence  or  mask 
the  aerodynamic  and  thermal  effects  being  tested.   As  the 
model  must  be  examined  after  the  test,  the  extent  of  radio- 
active contamination  as  a  biological  hazard  also  is  examined. 


Radiation  damage  is  caused  by  neutron  bombardment 

20 
and  interaction,  alpha  and  beta  bombardment  and  gamma 

radiation. 
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The  radiations  may  be  divided  into  two  groups.   The 
light  group  consisting  of  betas,  gammas  and  other  electro- 
magnetic radiations,  and  the  heavy  consisting  of  neutrons, 
fission  fragments,  accelerated  ions,  and  alphas.   As  a 

general  rule  the  damage  from  the  heavy  particles  is  more 

2  0 

severe  that  that  from  the  light  particles.     The  heavy 

particles  displace  the  atoms  of  a  solid  crystal  lattice 
from  their  normal  position.   In  general  this  will  result 
in  an  increased  ultimate  strength,  a  decrease  in  elongation 
and  a  reduction  in  area  of  the  material,  also  the  micro- 
hardness  of  the  material  increases  and  there  is  a  decrease 
in  impact  strength.   These  effects  are  minimized  if  the 
neutron  bombardment  occurs  at  elevated  temperature,  as 
is  the  case  in  this  test  facility. 

Graphite,  a  prime  candidate  for  a  re-entry  material, 
incurs  crystallite  damage  during  irradiation.   Hardness 
and  strength  are  increased  while  thermal  and  electrical 
conductivities  decrease.   Graphite  shapes  change  in  gross 
dimension  during  radiation.   All  of  these  effects  decrease 
in  magnitude  when  the  temperature  of  the  graphite  during 

radiation  is  increased,  with  almost  no  effect  when  graphite 

o 

temperatures  are  above  1000  C.   Also  the  integrated  flux 

levels  at  which  these  effects  begin  to  be  noticable  is 

19  20 

high,  of  the  order  of  10   neutrons  per  square  centimeter. 

As  this  test  facility  will  provide  very  high  model  tempera- 

IS 
tures  and  neutron  fluxes  at  the  model  of  the  order  of  10 

neutrons/  cm1  sec,  it  is  expected  that  irradiation  damage 

will  have  negligible  effect  on  the  test  results. 
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4.2  Biological  Hazards 

In  determining  the  biological  radiation  hazard  due 
to  model  contamination  the  following  assumptions  are 
made: 

i.   radioactive  particles,  such  as,  fuel,  fission 
fragments  do  not  adhere  to  the  model  upon 
impact, 

ii.  contamination  is  due  to  neutron-model 
interaction  causing  a  transmutation  to 
a  possibly  unstable  isotope, 

iii. contamination  is  acceptable  if  the  model 

can  be  handled  with  minimum  protection  within 
a  short  time  of  the  test  run, 


iv.  acceptable  radiation  levels  from  the  model 

_3 
are  7.5  x  10    rem/hour  as  recommended  by  the 

21 
US  National  Committee  on  Radiation  Protection. 


4.2a contamination  Environment 

Neutron  density  at  the  model  has  three  sources; 
direct  streaming  from  the  reactor  core,  neutron  flux 
from  the  fissionable  mass  in  the  test  section,  and  secondary 
emission  from  fission  fragments. 

Direct  streaming  of  neutrons  from  the  core  can  be 
measured  by  noting  the  attenuation  of  the  neutron  flux 
by  the  gases  in  the  nozzle  between  the  model  and  the  core. 
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The  differential  change  in  neutron  flux  due  to  absorption 
by  the  gas  is 


&  =  <?£.  -w.i) 


The  macroscopic  absorption  cross  section  of  the  gas 
varies  with  position  in  the  nozzle  as  the  density  of 
the  gas  decreases.   Assuming  this  variation  is  linear 
with  axial  nozzle  position,  %   ,  and  that  the  density 
goes  to  zero  at  the  model,  Equation   4.1   becomes 


&-*#£*.    (i-t)cl*  (4-2) 


When  integrated  and  evaluated  at  the  model,  the 
flux  at  the  model  is  given  by 

^   -  fa  e  (4.3) 

where  L  is  the  axial  position  of  the  model  and  Xo. 
is  the  absorption  mean  free  path 

The  axial  position  of  the  model,  \_j    yis  limited 
by  heat  transfer  considerations  in  the  nozzle  and  inlet 
and  must  be  minimized.   Then  L»   is  the  order  of  the 
required  nozzle  length  to  expand  the  flow  to  Mach  3, 
approximately  30  cm,  for  the  throat  size  set  in  Section  6.1 

For  the  chamber  condition  of  this  study  the  mean 
free  path  for  absorption  in  the  chamber  is  about  five 
times  this  value,  then  the  free  streaming  flux  from  the 
core  to  the  model  is  approximated  by 


$4  =  0.9  <&>  (4.4) 
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The  neutron  flux  from  the  mass  of  fissionable  material 
in  the  test  section  is  found  by  considering  the  change   in 
flux  as  the  flow  goes  from  the  critical  region  of  the 
core  to  the  subcritical  test  section. 

The  time  of  flight  from  the  core  to  the  test  section 
is  of  the  order  of  the  length  divided  by  the  speed  of 
sound  at  the  nozzle  throat, 


tf     -  U/C*  (4-5) 


If  the  test  section  is  considered  to  behave  as 
an  infinite  cylinder  with  constant  cross  section,  the 

critical  buckling  is  a  function  of  the  geometry  and  is 

9 

related  to  the  test  section  radius  by 

The  effective  multiplication  factor  in  the  test  region 
is  i 


>  +(^r  h.       <4-7) 

where  ^^   is  defined  by  Equation  3.9 

With  the  assumption  that  the  number  of  neutron  generations 
is  given  by  the  time  of  flight,  *■$.  f      divided  by  the  average 
neutron  lifetime,  X  a./  il*  ,  the  neutron  flux  at  the  model,  due 
to  the  reactivity  of  the  fissionable  gas  in  the  test  section, 
is  given  by 

_tt(lu«-ifc  (4.8) 
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Equations    4.5    through    4.3    are  evaluated  for 
the  conditions  in  this  facility, 

ft 

keg'   .002 
5©C 


_  4 

-ts  -   5*  /D       seconds 


fafcJL  ». 


to  obtain  the  flux  at  the  model  due  to  the  reactivity  of 
the  test  section,  -ISO 

^k  -  ii  e  (4.9) 

As   ^yic     is  very  much  smaller  than  <rf*  it 

is  considered  negligible  in  this  analysis. 

Secondary  emission  of  neutrons  by  fission  fragments 
in  the  vicinity  of  the  model  will  also  contribute  to  the 
neutrons  bombardment  of  the  model.   Given  a  fission  fragment, 
A  ,   with  a  half  life,  "ty^A  /   and  a  yield  fraction 
per  fission  of  /?A  #  from  the  definition  of  half  life,  the 
radioactive  decay  constant  #  A  ,  is 

>A  -   J"12-  (4.10) 

The  rate  of  creation  of  atoms  of  fragment   A   is 

iA  --  "%Z%4o  (4.11) 

dt 
where  £,^9£ is  the  number  of  fissions  per  unit  volume  per 

unit  time. 

The  rate  of  decay  is 

d*   -  -Xa  k  (4.12) 

Thus  in  the  reactor,  fission  products,  A  ,  build  up  at 


the  rate 


i|  =  *Mi,&  "***  (4.13) 
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and  concentration  of  species  A    at  time  t  is 

Aftl  -  £  (% ^i4  -  fot^. -»k] e^)  (4.14) 

With  the  initial  condition  that  the  initial  concentration 

of  species  A  is  zero,  the  concentration  of  fragments  leaving 

the  core  is  given  by  A(to')  ,  where  to   is  the  dwell  time 

of  fluid  in  the  core, 

{.     *  ?»^?  (4.15) 

At   t&  =lseco»\J,  the  concentration  of  A  leaving  the 
reactor  is  given  by 

A(fl.  n»£*0o  (i-e**)  (4.16) 

The  rate  of  decay  in  the  nozzle  is  given  by  Equation  4.12 
evaluating  this  at  the  test  section  with  a  time  of  flight, 
t$  ,     from  the  core,  the  rate  of  decay  of  fragment  A 
at  the  model  is 


J*  *  -  7u  ± *»  *•  ( i_  eA*)  e~x**f  <4-17> 


If  it  is  assumed  that  one  neutron  is  emitted  per 
decay,  the  sum  over  all  fission  fragments  of  the  decay 
rates  given  by  Equation   4.17  is  the  rate  of  creation 
of  neutrons  at  the  model, 

^-  2>e**(i-e**^.j*.  (4.i8, 

A 

The  absorption  mean  free  path  in  the  test  section 
is  about  fifty  times  larger  than  the  test  section  radius, 
therefore  all  neutrons  freed  in  the  test  section  will  hit 
the  model  or  walls.   If  all  neutrons  hitting  the  walls 
are  absorbed,  the  ratio   of  neutrons  hitting  model  to  those 
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hitting  the  wall  is  given  by  the  ratio  of  model  radius  to 
test  section  radius,   Rw  /  Rt       .   For  negligible 
model  blockage  of  the  test  section  assume  this  ratio  is  0.1 

Then  the  delayed  neutrons  strike    a  unit  area  of 
the  model  at  a  rate 


21TR-  R*   cU  (4.19) 

The  fission  products  yield  are  given  in  Reference  22. 

An  examination  of  the  distribution  of  fission  products 
shows  that  for  a  mass  number  less  than  80  or  greater  than 
160,  Y{a  is  less  than  one  part  in  ten  thousand  and  the 
contribution  to   ^S^   is  negligible  regardless  of  the  half 
life.   Also,  if  "ti/x  is  greater  than  one  minute,  XA  is  less 
than  .01  seconds   ,  this  combined  with  a  maximum   Y7*   of 
ten  percent  from  Reference  2  2 ,      shows  that  the  contribution 
of  any  element  with  a  half  life  much  greater  than  one  minute 
may  be  ignored  in  approximating  the  neutron  flux  due  to 
secondary  emission  from  fission  fragments.   References 20   and 

2  2   were  examined  for  fission  products  with  yield  fractions 

-4 
greater  than  10    and  half  lives  less  than  one  minute.   The 

isotopes  satisfying  these  criteria  are  given  in  Table  1.   With 

these  values,  the  rate  of  creation  of  neutrons  at  the  model  is, 

-2 

by  Equation  4.18,  1.72  x  10    neutrons  per  second  and  the 

neutron  flux  at  the  model  due  to  secondary  emission  of  neutrons 
is 

**""»»«**.  (4.20) 

It  is  noted  that  this  method  gives  a  conservatively 
high  estimate  of  flux  at  the  model  due  to  assuming  that 
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the  total  fission  yield  for  a  given  atomic  number  indicated 
in  Reference  23  is  the  element  with  the  atomic  number  with  a 
short  half-life,  and  further  that  these  radionuclides  all 
emit  a  thermal  neutron  upon  decay. 

There  are  two  unconservative  assumptions  in  this 
development,  that  is  that  the  initial  condition  is  the 
absence  of  fission  fragments  at  the  beginning  of  the  reactor 
dwell  time,  and  that  only  short  lived  neutron  emitters  are 
involved.    In  operating   as  a  closed  cycle  , the  supply  of 
nitrogen  and  uranium  will  be  contaminated  by  the  fission 
fragments  generated  during  previous  residence  in  the  reactor 
core,  increasing  the  initial  concentration  of  fragments  and 
allowing  build  up  of  fission  fragments  with  a  long  half-life. 

Due  to  the  shear  magnitude  of  the  power  involved  in 
this  facility,  it  is  expected  that  the  time  for  a  given  mass 
of  fluid  to  complete  a  cycle  will  be  long  enough  to  allow 
decay  in  the  concentration  of  the  prompt  neutron  emitters 
and  that  provision  can  be  made  for  monitoring  and  removal 
of  long  lifed  fission  fragments  as  a  part  of  the  fuel  handling 
process . 

With  these  assumptions,  the  flux  due  to  secondary 
emission  from  fission  fragments  is  several  orders  of  magnitude 
smaller  than  the  direct  streaming  flux,  and  the  former  is 
neglected. 

The  flux  striking  the  model  is  then  approximated  by 
the  direct  streaming  neutron  flux  of  Equation  4.4  . 
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The  average  flux  in  the  gaseous  core  is  set  by  the 
reactor  size,  the  power  required,  the  energy  per  fission,  €, 
and  the  macroscopic  cross  section  for  fission  of  the  fuel  in 
the  reactor  core, 

0o   -  YT~  (4*21) 

At  a  power  level  of  7100  megawatts  with  a  1.7  m  core, 

15         ,   2 
the  average  flux  in  the  gas  core  is  3.26  x  10    neutron/cm   sec 


Equation  4.4  gives  the  flux  at  the  model -for  these 

.  .  15  2 

conditions,  2.93  x  10    neutrons/  cm   sec. 


4.2b  Model  Contamination 

The  radiation  hazard  from  the  model  is  assumed  to  consist 
of  secondary  emission  from  neutron  flux  induced  radionuclides 
on  the  model.   The  governing  equation  is  of  the  same  form 
as  the  build  up  of  radionuclides  in  the  core,  then  decaying 
after  exposure.   The  build  up  of  species  A  on  the  model  during 
a  test  run  of  time,  t-t   ,  is 

a(M  -  ^X0"  ('-e^')  (4-22) 

where  <£•»,*   is  the  neutron  macroscopic  absorption  cross  section 

in  the  reaction  b+o*"*1  — *A  and  B  is  one  of  the  original 

materials  of  the  model  or  a  transmuted  species.   The  rate  of 
decay  of  A  after  the  test  is 

^-^tfuO-^^e***  <4-23) 

where  t  is  the  time  since  exposure  to  the  neutron.   After 
the  decay  rate  for  a  species  has  been  determined,  the  type 


38 

of  decay  must  be  examined  to  determine  the  mode  of  decay,  that 
is  alpha,  beta,  gamma  or  neutron  emission  and  the  energy  of  the 
emitted  particle.   The  particle   type  and  energy  set  the  maximum 
decay  rate  that  a  man  can  be  exposed  to  and  not  exceed  the 
limits  of  .006  rem/day.   It  is  noted  that  if  more  than  one 
species  is  formed  or  more  than  one  mode  of  decay  is  present, 
the  sum  of  the  exposures  will  determine  the  allowable  decay 
rates  of  each  species  or  mode. 


As  graphite  is  the  prime  candidate  for  material  to  be 
tested,  it  will  be  examined.  Carbon  undergoes  the  reactions 


*c* 

+ 

1 

•n     = 

(,C 

stable    • 

(4.24a) 

*c" 

+ 

on'     = 

.c- 

ti/x    =5570   years 

(4.24b) 

tc" 

+ 

•n1     = 

/3    +   tn 

stable 

(4.24c) 

22         3 
for  KlclX   =  8.55  x  10    atom/cm 

r-      ,  -   ,rt-27   2 
veu  c,i  =  3.3x10     cm 

-27    2 
e*mm       =  .9  x  10    cm 

tT    =  200  seconds 

The  initial  decay  rate  for  the  model  is,  by  Equation  4.23 

With  no  absorption  in  the  model,  the  surface  beta  flux 
for  a  cylinderical  model  is  conservatively  given  by 

dC„  Vow.,  S  clCu   g^ 

dt   Am      oit    2.  (4.26) 

For   a  carbon  model  with   a   radius    of    1.25   centimeters 


d"t  tfma   Sec 
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Reference  20  gives  the  exposure  of  betas  of  various 
energies  which  produces  an  allowable  dose  of  .006  roentegen 
equivalent     man  per  eight  hour  day.   The  exposure  rate 
above  is  on  the  safe  side,  therefore  a  graphite  model  may 
be  handled  immediately  after  a  run  at    7,100      Megawatts 
for  200  seconds. 


Correlations  between  nuclear  decay  products  and  allowable 

9,20,23 

dose  rates  are  readily  available  in  the  literature 

A  very  slight  protective  layer,  such  as  paper  or  cloth,  will 

provide  protection  against  alpha  particles,  which  do    not 

9,20 

present  an  external  radiation  hazard. 


It  is  also  probable  that  high  percentages  of  hydrogen  2 
will  be  present  in  an  ablative  material,  and  the  decay  of 
hydrogen  isotopes  yields  a  low  energy  beta  particle  in 
concentrations,  determined  by  Equation   4.22  l    five  orders  of 
magnitude  below  carbon  and  therefore, as  a  graphite  model, 
also  may  be  handled  immediately  after  a  test. 

This  analysis  can  not  possibly  determine  if  each  and 
every  material  will  be  safe  to  handle  after  a  test.   The 
user  will  have  to  examine  the  reaction  of  each  model 
individually  to  determine  the  contamination  from  each  material 
and  the  time  history  of  resulting  total  model  contamination, 
as  given  by  Equation  4.22  summed  over  all  model  components. 
The  resulting  model  activity  can  be  compared  to  allowable 
dose  rates  to  determine  if  and  when  the  model  can  be  handled. 

It  is  noted  by  examination  of  Equations   4.4,  and  4.22 
that  the  decay  rate  at  the  model  is  directly  proportional  to 
the  flux  rate  in  the  reactor.   This  in  turn,  by  Equation  4.20, 
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is  dependent  on  the  specific  power  of  the  reactor  and  varies 
inversely  with  the  concentration  of  uranium  in  the  core,  so 
that  low  power  density  and  high  uranium  concentration  are 
desired . 

As  is  shown  in  Section  6.1,  the  operating  power  of 
the  reactor  is  proportional  to  the  area  of  the  core,  with 
core  temperature  and  pressure  fixed.   The  specific  power  varies 
then  inversely  as  the  core  radius,  as  does  the  concentration 
of  uranium  in  the  core  at  small  radii,  shown  in  Figure  7. 

It  can  be  concluded  that,  at  least,  for  many  possible 
re-entry  materials  there  would  be  no  difficulty  in  handling 
the  models  after  testing  and  that  model  contamination  is  not 
a  critical  parameter  in  the  design  of  this  system,  but  should 
a  particular  material  to  be  tested  become  highly  contaminated, 
the  degree  of  contamination  can  be  reduced  by  reducing  either 
the  test  duration  or  the  specific  power  density.   It  is  assumed 
that  this  does  not  occur  in  this  study. 


41 


CHAPTER  V 

REACTOR  COOLING 

5.1  Heat  Loads 

The  walls  of  the 'gaseous  core  reactor  are  subjected  to 
high  heat  loads  due  to  convective,  conductive  and  radiation 
heat  transfer  and  heat  producing  nuclear  reactions  in  the 
solid  materials  of  the  core  materials. 

At  the  high  temperatures  of  the  core  the  heat  load  due 
to  radiation  transfer  is  much  greater  than  the  convective 
and  conductive  heat  loads:  the  convective  and  conductive 


heat  loads  are  considered  negligible  in  this  analysis.   The 
core  is  assumed  to  radiate  as  a  b! 
the  walls  absorptivity  set  as  0.6 


o 
core  is  assumed  to  radiate  as  a  black  body  at  10,000  K  and 


The  heat  load  due  to  nuclear  reactions  is  assumed  to 

be  due  to  neutron  and  gamma  fluxes  from  the  core  and  is  set  as 

9 

eight  percent  of  the  power  output  divided  by  the  core  area- 


The  heat  load  per  unit  wall  area  is 

Qoc-  6  <rTo4*  .08  P/Ac  (5-1) 

where 

6  is  the  wall  absorptivity,  and 

(J*  is  the  Stephan-Boltzmann  constant 
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For  the  conditions  of  this  study,  the  heat  load  to  the 
wall,  given  by  Equation  5.1   is 


0^=  S.4B  >K>3  **ty 


Cm  -  5€C. 


5.2  Criticality  Effects 

The  loads  given  by  Equation  5.1   are  much  too  high  for 
the  thick  beryllium  moderator  to  conduct  away  from  the  core, 
and  are  predominantly  due  to  radiation.   Reference  7  suggests 
the  use  of  a  thermal  barrier,  cooled  by  transpiration 
and  radiation  seeding  to  protect  the  walls  of  a  gaseous  core 
nuclear  rocket  under  similar  high  radiation  heat  loads. 
This  method  of   cooling   is  not  appropriate  for  use  in  the 
reactor  used  in  this  system.   The  large  size  of  the  reactor 
is  due  to  the  high  absorption  cross  section  of  the  working 
gas, nitrogen.   This  causes  a  low  value  of  the  infinite 
multiplication  factor  in  the  reactor,   Koo  ,  given  by 
Equation   3.9  .   Examination  of  Equation  3.9  shows  that 
increasing  the  macroscopic  absorption  cross  section  in  the 
core  will  reduce  the  infinite  multiplication  factor,  K.»  . 
As  this  is  already  very  close  to  1.0  in  this  reactor,  addition 
of  even  a  small  thermal  barrier  causes  the  reactor  to  be 
subcritical  for  any  reactor  size. 

Even  if  a  material,  with  sufficiently  low  neutron 
absorption  characteristics  and  high  strength  could  be  found, 
the  transpiration  flow  necessary  to  internally  cool  the 
thermal  barrier  is  not  available  due  to  the  large  area  of 
this  core  and  the  relatively  small  throat  area. 
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The  moderator  in  this  study  is  protected  from  the 
high  thermal  radiation  heat  loads  by  a  seeded  layer  of  carbon 
particles  in  the  core  near  the  walls.    The  small  amount  of 
'carbon  necessary  to  protect  the  walls   does  not  preceptibly 
lower  the  multiplication  factor,  but  has  a  more  pronounced 
effect  on  lowering  the  diffusion  length,  defined  by  Equation 
3.8    ,  so  that  the  carbon  acts  as  a  moderator  and  slightly 
decreases  the  critical  size,  as  shown  by  Equation   3.10 
This  effect  is  due  to  the  very  small  absorption  cross  section 
of  carbon,  which  effects  K»  ,  and  the  larger  scattering 
cross  section,  which  reduces  L»   .    The  effect  of  the  carbon 
seed  on  criticality  is  ignored  in  this  analysis. 

5.3  Reactor  Flow  Pattern 

In  order  to  protect  the  moderator  of  the  core,  it  is 
necessary  to  establish  the  following  flow  pattern.  Consider 
any  diagonal  plane  which  includes  the  axis  of  symmetry  of 
the  nozzle.  The  reactor  core  flow  and  the  nozzle  are  axi~ 
symmetric.  The  reference  direction, 9  -  O ,  is  aligned  with 
the  nozzle  axis;  the  flow  has  a  source  at  0'TT  ,  the  fuel 
and  nitrogen  inlet,  and  a  sink  at  the  nozzle.  A  boundary 
layer  is  formed  along  the  walls  as  the  flow  proceeds  from 
inlet  to  the  nozzle,  much  like  a  low  speed  divergent-con- 
vergent nozzle. 

The  uranium  fuel, in  the  form  of  a  fine  powder  mixed 
with  nitrogen,  is  sprayed  into  the  core  at  the  0  -=•  ^f 
position.   The  nitrogen  is  sprayed  into  the  core  through  a 
annulus  around  the  fuel  nozzle.   Perfect  mixing  of  the  fuel 
and  nitrogen  is  assumed.   It  is  assumed  that  the  carbon  seed 
tends  to  remain  toward  the  walls  and  is  limited  to  a  layer 
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one-tenth  of  the  core  radius  thick  and  starting  at  the  core 
wall. 

The  gamma  and  neutron  fluxes  to  the  moderator  is  un- 
affected by  the  thermal  radiation  barrier.   The  total  heat 
load  to  the  moderator  is  the  sum   of  the  nuclear  heat  load 
and  the  thermal  radiation  which  is  not  absorbed  by  the  seeded 
layer.   The  moderator  is  cooled  by  a  combination  of  radial 
injection  of  a  fraction  of  nitrogen  flow  into  the  core,  and 
water  cooling.   The  nitrogen  flow  itself  is  not  sufficient 
to  cool  the  moderator  and  must  also  be  used  to  establish  the 
desired  flow  pattern. 

5.4  Radiation  Barrier 

With  the  flow  pattern  in  the  core  as  established  in 
paragraph  5.3,  the  flow  is  essentially  flow  in  a  nozzle 
with  a  diverging  section  proceeding  the  converging-diverging 
nozzle.   Reference  12  gives  the  results  of  using  a  seeded 
layer  to  protect  a  gaseous  core  rocket  nozzle.   With  perfect 
mixing  of  the  fuel  and  working  fluid,  a  seeded  layer  1/10 
of  the  local  nozzle  radius  thick  will  reduce  the  radiant 
energy  flux  by  two  orders  of  magnitude  if  the  seeded  layer 
absorption  coefficient,   K.   ,  is  1/3  centimeter. 

Reference  13  gives  a  correlation  of  the  total  absorption 
coefficient  per  seed  particle  concentration,  K  /** «  and  the 
particle  radius.   Using  carbon  particles  as  the  seeding 
material,  it  is  desired  that  the  average  concentration  of 
carbon  particles  be  minimized,  both  for  nuclear  and  thermo- 
dynamic reasons. 
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The  average  carbon  density,  0   ,  in  the  seeded  layer  is 


?    -fetUftf) 


where 

ft      is  the  mass  density  of  the  carbon  seed 
and    R.a     is  the  average  seed  radius 

,,   -» 

With  K.  set  at  As   Cm  /  the  average  density  is  minimized 

at  3.T  *  10   c\ ram /c«ntL meter  *,  this   occurs  with  a  concentration 

\o  3 

of  4.5  k  io       particles/cm  having  an  average  diameter  of 

0.1  microns. 


As  the  seeded  layer  does  not  fill  the  core,  the  average 
seed  density  over  the  core  is 


f-jff(R?-9R?)  ji.3, 

Je  4 


|irRc3 

=-.    j_tflO      <ara.vn4/ce«ti.me.^er 

The  contribution  to  the  macroscopic  absorption  cross 
section  due  to  the  carbon  atoms  in  the  core  is  given  by 

where  CA.      is  Avogadros  number  and c/\c is  the  molecular  weight 
of  carbon* 

For  the  concentration  of  carbon  in  the  core  due  to  the 


seeded  thermal  barrier,  Equation  5.4yields  an  absorption 
cross  section  contribution  of  \.G  *  \0   ce-i.U»-/ieVer~      .   i 
this  is  four  orders  of  magnitude  less  than  the  absorption 
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cross  section  of  the  nitrogen  and  uranium  mixture  in  the 
core,  the  presence  of  the  -carbon  seed  on  criticality  is 
ignored . 

5.5  Moderator  Cooling 

The  seeded  radiation  barrier  discussed  in  Section  5.4 
reduces  the  radiant  heat  load  to  the  wall,  given  in  Section 
5.1,  by  two  orders  of  magnitude,  but  does  not  effect  the  nu- 
clear radiation  induced  heat  load.   The  resulting  heat  load 

2 

to  the  moderator  is  350   cal/cm   sec. 

The  thickness  of  the  moderator  and  the  thermal  conductivity 
of  beryllium  prevent  this  heat  load  from  being  conducted 
through  the  moderator  and  removed  by  back  side  cooling. 

A  portion  of  the  incoming  nitrogen  flow  is  radially 
passed  through  the  moderator  to  absorb  the  heat  load  to  the 
moderator.   The  coolant  nitrogen  flow  is  injected  tangentially 
into  the  core  in  such  a  manner  as  to  be  parallel  to  the  flow 
of  the  boundary  layer  at  the  point  of  injection.   This  is 
to  prevent  the  injected  flow  from  driving  the  seeded  thermal 
barrier  away  from  the  wall.   A  small  amount  of  carbon,  as 
is  necessary  to  replace  carbon  in  the  thermal  barrier  which 
vaporizes  or  diffuses  out  of  the  thermal  barrier,  is  carried 
by  the  cooling  nitrogen  flow. 

The  maximum  temperature  of  the  beryllium  moderator  is 
set  at  1000  K  by  structural  strength  limits -With  one  half 
of  the  primary  nitrogen  flow  used  to  establish  the  desired 
flow  pattern,  the  remaining  primary  nitrogen  flow  is  capable 
of  absorbing  one  third  of  heat  load  in  the  moderator,  increasing 
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in  temperature  from  293  K  to  1000  K.    The  remaining  heat 
load  is  absorbed  by  bleeding  four  percent  of  the  water  flow 

enroute  to  the  water  diffusor  and  allowing  its  temperature 

o 

to  increase  to  1000  K  in  the  moderator. 
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CHAPTER  VI 
NOZZLE  HEAT  TRANSFER 

6.1  Nozzle  Size 

The  hot  gases  in  the  core  are  radiating  energy  away 
from  the  core.   As  the  moderator  walls  are  not  capable  of 
transferring  this  heat  load,  they  must  be  protected  by  a 
seeded  thermal  barrier ,  and  the  absorbed  energy  carried  out 
of  the  core  through  the  nozzle.   Thus  the  minimum  operating 
power  is  set  by  the  radiant  heat  load  from  the  core,  and  the 
size  of  the  core  wall  surface.   The  nozzle  throat  area  is 
found  by  equating  the  radiative  power  in  the  core  to  the 
energy  flow  through  the  nozzle,  then, 

A*  ■  TOT 

For  the  operating  conditions  set  in  Section  2.2,  the 

2 

minimum  throat  size  is  93  centimeters  . 

6.2  Heat  Loads 

Large  heat  loads  are  experienced  by  the  walls  of  the 
nozzle  as  the  flow  is  accelerated  from  stagnation  conditions 
of  the  core  to  the  high  Mach  number  desired  in  the  test 
section.   The  heat  load  is  due  to  convective  and  radiative 
heat  transfer  from  the  hot  gases  and  nuclear  heating  of  the 
nozzle  walls  by  alpha,  beta  and  gamma  interaction.   To 
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prevent  failure  of  the  nozzle,  this  heat  load  must  be 
transferred  through,  absorbed  by  or  blocked  from  the  nozzle 
walls. 

Turbulent  flow  is  assumed  in  the  nozzle,  also  it  is 
assumed  that  the  flow  is  frozen, isentropic  and  one  dimensional 
in  determining  local  flow  properties.   The  heat  generated  in 
the  nozzle  by  nuclear  interaction  is  assumed  negligible  compared 
to  the  convective  and  radiation  heat  loads. 

6.3,  Nozzle  Configuration 

The  nozzle  consists  of  a  series  of  longitudinal  porous 
tubes,  the  area  ratio  of  which  is  adjusted  to  provide  an 
over  pressure,  of  the  cooling  fluid  in  the  nozzle  tubes 
above  the  pressure  of  the  flow  in  the  nozzle,  of  10  atmospheres. 
This  provides  the  driving  potential  for  transpiration  cooling 
of  the  nozzle.   Flow  inside  the  coolant  tubes  is  assumed  to 
be  turbulent. 

The  tube  cross  section  is  a  "U"  shape,  the  flat  sides 
of  which  may  vary  to  provide  the  desired  area  ratio  and  also  to 
facilitate  joining  the  tubes  together  to  form  a  smooth  nozzle 
shape.   The  heat  transfer  surface  is  approximated  by  a  circle 
whose  radius  is  R.   The  fin  effect  of  the  radial  sides  of 
the  cooling  tubes  is  not  accounted  for,  giving  a  conservative 
value  for  the  heat  transfer  to  the  coolant. 

The  coolant  tubes  support  the  thermal  stress  due  to 
heat  transfer  across  their  radially  inner-most  face.   They 
also  support  the  hoop  stress  due  to  the  coolant  pressure. 
The  tubes  are  closed  by  a  shell  which  supports  the  pressure 
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of  the  hot  gas  in  the  nozzle. 

6.4  Maximum  Heat  Transfer 

The  maximum  heat  load  that  can  be  transferred  through 
the  wall  is  set  by  the  material  properties,  the  configuration 
of  the  wall  and  the  pressure  the  wall  must  support. 

For  a  thin  walled  tube,  with  the  wall  thickness,  d,  very 
much  less  than  the  tube  radius,  R^ ,  the  temperature  gradient 
across  the  tube  is  linear.   The  thermal  stress  across  the 
tube,  in  the  tangential  direction,  are  maximum  in  tension  at 
the  cool  inner  fiber, 

<h*  "  *aEL0  (Tw-T-o)/e  (6-2a) 

and  maximum   in   compression   at    the  hot    outer    fiber 

<Hv>  =     -*EL      (Tv-TwoVZ.  (6.2b) 

where 

CC  is  the  coefficient  of  linear  expansion, 

£  is  the  modulus  of  elasticity, 

lit   is  the  hot  side  wall  temperature, 

"T^is  the  cold  side  wall  temperature, 

and  both  0C  and  E  are  functions  of  temperature. 

The  pressure  differential  across  the  cooling  tube  wall, 
-p,  causes  a  uniform  hoop  stress  in  the  wall, 

a  (6.3) 

where 

d.  is  the  wall  thickness, 
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and  the  positive  sign  indicates   tensile  hoop  stress,  the 
negative  compressive. 

The  total  stress  is  given  by  the  sum  of  the  thermal 
and  hoop  stress;  scaler  sum  in  this  case  as  both  are  tangential 

At  the  extreme  fibers, 

CT=r  1  QCECT.-Two]     +     1°Rt  (64) 

2  d 

The  configuration  chosen  in  this  analysis  utilizes  a  series 
of  small  cooling  tubes  forming  the  nozzle.   This  gives  a 
tensile  hoop  stress  and  the  maximum  total  stress  is  tensile  at 
the  cool  inner  fiber.   The  hoop  stress  relieves  the  compressive 
thermal  stress  at  the  hotter  inner  wall,  where  the  material  is 
weakest. 


The  heat  transferred  through  the  wall  is  given  by, 

^  =  df      ~      '  (6.5) 


where 

k  is  the  wall  thermal  conductivity 

Equation  6.4  is  solved  for  the  temperature  across  the 
wall,  and  this  is  substituted  into  Equation   6.5   , 

0  =      -    f (6.6) 

Zh 

The  maximum  heat  transfer  through  the  wall,  with  respect  to 
wall  thickness, d,  is  obtained  by  setting,  c\Q/dcL  equal  zero 

«Ed's  *edl"    u  <6-7> 
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The  condition  for  maximum  heat  transfer  through  the 
wall  is, 

.    ^  .  -e«*  -     .e;  -   «„  (6.8) 

The  maximum  heat  flow  through  the  wall  is 

Q-MM  *     2^fR,  (6*9) 

Here  failure  bccurs  when  the  tensile  stress  at  the  inner 
fiber  of  the  coolant  tubes  reaches  the  ultimate  stress  of  the 
tube  material.   Radial  stress  is  due  to  the  pressure  differen- 
tial across  the  tube  thickness  and  is  negligible  compared  to  the 
tangential  stress.   The  nozzle  has  incroporated  an  expansion- 
slip  joint  at  the  up  stream,  low  Mach  end;  this  prevents 
longitudinal  thermal  stress  so  that  the  stress  in  the  coolant 
tubes  is  considered  uniaxial. 

6.5  Material  Selection 


a 


A  material  was  sought  with  a  maximum  value  of  -j"^ 

in  order  to  maximize  heat  transfer  with  respect  to  material 
selection.   This  parameter  is  a  function  of  temperature  in  a 
given  material  and  the  temperature  must  be  high  enough  so 
that  super-cooling  is  not  necessary  at  the  back  side  to  absorb 

the  heat  load.   Wall  temperature  also  affects  the  heat  load 

10 

to  the  wall  but  examination  of  the  Bartz's  Equation    , 

utilized  to  determine  the  convective  heat  load,  shows  only 

a  weak  dependence  of  convective  heat  transfer  on  wall  temperature 

24  o 

This  investigation  chooses  Rene  41    at  Tw  of  1025  K 

as  having  the  desired  structural  properties. 
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At   this   temperature, 

2 

GtL  is    146,000'  pounds/inch 

—  7  2 

h     is    2.4   x    10      pounds/inch 

•  o 

K  is  0.05  cal/cm-sec-  K 

Ok  1.57  x  10    cm/cm  K 

Minimum  wall  thickness  is  set  by  manufacturing  limitations 
and  is  assumed  to  be  0.01  cm  for  this  analysis. 

This  choice  of  material  and  wall  thickness  gives  a 
maximum  heat  transfer  per  unit  wall  area  of 

Qm„  »  1.94  *»0  <*aX/cm*  sec 

provided  that  the  cooling  tube  radius,  set  by  Equation 
6.8  with  a  ten  atmosphere  pressure  differential  across  the 
wall,  is  5  cm  or  less.   The  temperature  drop  through  the 
wall  is 

T~  -  Tw.   -  ^89%. 

and  TWo       =  G3C*k. 

6.6  Test  Section  Mach  Number 

The  nozzle  shape  considered  is  a  convergent-divergent 
nozzle  formed  by  the  previously  defined  cooling  tubes. 
The  nozzle  is  axially  symmetric  with  a  45  degree  convergent 
section,  15  degree  divergent  section  and  a  radius  of  curvature 
at  the  throat  equal  to  the  throat  inner  radius. 

The  heat  load  to  the  nozzle  is  weakly  dependent  upon 
throat  radius,  affecting  the  convecting  heat  transfer  co- 
efficient only  in  this  analysis.   The  pump  work  required  of 
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the  system  as  a  whole ,  however,  varies  directly  with  the 
mass  flow,  which  varies  as  the  area  of  the  throat. 

The  throat  radius  is  chosen  in  Section  6.1  as  5.45 
centimeters . 

The  test  section  Mach  number  is  arbitrary  with  a  friction- 
less  isentropic  nozzle.   In  this  nozzle,  transpiration  cooling 
is  required  until  the  flow  has  expanded  to  about  Mach  6, 
which  requires  a  lengthy  nozzle,  of  the  order  at  15  meters. 
Expanding  the  flow  to  Mach  6  increases  the  transpiration 
coolant  requirement  and  nozzle  losses.   For  this  analysis,  a 
test  section  Mach  number  of  3.0  is  chosen  to  enable  simulation 
of  the  desired  flow  fields,  but  reduce  transpiration  pumping 
requirements  and  nozzle  losses. 

6.7  Transpiration  Cooling  Requirements 

The  heat  load  per  unit  area  to  the  wall,  neglecting 
gamma  heating,  is  given  by, 

Q   "  Oo„  *Qr  (6.10) 

where  Q©<.  is  the  convective  heat  load  in  the  absence  of 
transpiration  cooling, 


Qr  is  the  radiative  heat  load,  assumed  unaffected 
by  transpiration  cooling. 

The  convective  heat  load  was  obtained  by  the  method  of  Bartz 

for  turbulent  compressible  flow, 

Q-c-VxtT^-X)  t5-11*' 

where .   r  ^ A  ,  /  ,,o^  \       t    -  _  \,%   /  in  \ .  (~\  I  A  \  -  7 


10 


wi-m  (^mri  &) 


6"   (6.11b) 
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cr=   1 

[k    %(l*^')^]-8"»  [|*^rf]*    (6.11c) 

K  is  the  heat  transfer  coefficient, 

(J  is  0.6, 

D  is  the  throat  diameter,  inches/ 

Jl  is  the  viscosity, 

Cp  is  the  specific  heat  at  constant  pressure, 

"T"  is  the  temperature,  °R     , 

sub         =    wall, 
w 

sub         =    bulk, 
o  ' 

sub         =     adiabatic  wall, 
aw 

-f£  is  the  chamber  pressure, 

Yi    is  the  radius  of  curvature  of  throat,  inches, 
1^  is  1  for  the  nozzle, 
Ay£\  is  the  nozzle  area  ratio, 
y  is  the  ratio  of  specific  heats, 
M  is  the  Mach  number. 


14 
The  thermodynamic  properties  of  nitrogen   were  used  in 

this  analysis,  and  the  presence  of  the  uranium  was  ignored. 


The  radiation  heat  load  is  obtained  by  considering  the 
gas  to  radiate  at  stream  ambient  temperature  as  a  black  body 
and  the  wall  to  have  an  absorptivity  of  0.6, 

Q**€w<rTi  -~*w£.-g*~0  (6.12) 

The  radiation  from  the  wall  is  negligible  compared  to  the 
load  to  the  wall  and  this  reduces  to 


Qk   -  £w  6"  TA  (6.13) 
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where 

£vo  is  the  wall  absorptivity^ 

£„  is  the  gas  absorptivity, 

<T*  is  the  Stef an-Boltzmann  constant, 

T   is  the  Temperature,  ambient  free  stream  *k.y 

&w  is  the  wall  emissivity. 

The  heat  load  to  the  nozzle  wall  in  the  absence  of 
transpiration  cooling  is  given  in  Figure  8     .   The  maximum 
heat  load  that  the  wall  is  capable  of  withstanding  is  lower 
than  the  load  experienced.   Transpiration  cooling  by  the 
injection  of  low  temperature  nitrogen  is  used  to  reduce  the 
convective  heat  load  to  the  wall. 

The  effects  of  transpiration  cooling  in  reducing  the 
convective  component  of  heat  transfer  is  given  in  Referencell. 
The  results  are  given  in  terms  of  the  Stanton  number  reduction 
obtained  by  injection  of  a  homogeneous  coolant  into  a  turbulent 
nozzle  flow  verses  a  blowing  parameter,  F/st0.  Here  F  is  de- 
fined as  the  ratio  of  the  transpirant  mass  velocity  to  the 
local  free  stream  mass  velocity, 

V*  =  -pv"  (6.14) 

where  G  is  defined  as  the  transpiration  mass  velocity.  The 
Stanton  number  is  defined  as, 

U  (6.15) 

Then  the  Stanton  number  reduction  is,  to  the  first  order, 
the  reduction  in  convective  heat  load, 


St/St0  as  Q.c/QCo  .(6.16) 
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The  maximum  heat  load  transferred  through  the  wall, 
with  transpiration  cooling,  is  the  sum  of  the  reduced  convective 
heat  load  and  the  radiation  heat  load  minus  the  heat  absorbed 
by  the  coolant  passing  through  the  wall. 

CL..  «  Qc+Q*-^ep(Tw-Tw.,>       <6-17> 

where   Co         is  evaluated  at  local  pressure  and  the  average 
temperature  in  the  wall. 

The  convective  heat  load  without  transpiration  cooling 
is  given  by, 

0c0  =  Vw(TL.-T.)  <6-18) 

and 

¥*•-£-  ¥■     "  '   (6-19) 

With  Equations  6.18  and  6.19  the  maximum  heat  flow  through 
the  wall  can  be  rearranged  to 

The  left  hand  side  of  this  equation  is  known  at  each 
station  and  the  right  hand  side  suggests  the  possibility  of 
a  graphical  solution  using  the  difference  between  the  Stanton 
number  reduction  data  of  Reference  11  and  a  family  of  straight 
lines  of  slope  Cp  (TV--7N-*')  /CP  (Taw  -Tw  )         .   Figure  9 
reproduces  the  Stanton  number  reduction  data  of  Reference  11 
and  illustrates  the  graphic  method  of  solution.   The  blowing 
parameter,  F,  is  obtained  by  this  method  and  from  It,  the 
transpiration  mass  velocity  necessary  to  reduce  the  heat 
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load  through  the  wall  to  the  maximum  allowable.   This  flow 
is  graphically  integrated  to  obtain  the  total  required 
transpiration  flow. 

The  transpiration  cooling  flow  per  unit  axial  nozzle 
length,  2.TCRGi  /  is  given  in  Fi  gure  10   ,  along  with  the  total 
transpiration  mass  flow. 

Under  the  conditions  of  this  study,  the  transpiration 
coolant  flow  necessary  to  prevent  thermal  damage  to  the 
nozzle  and  test  section  is  I.Z5x  lO  Ara.m/secondt    ,  and 
the  ratio  of  transpiration  flow  to  the  primary  flow  in  the 
nozzle  is, 

a*  =  0.37 

This  estimate  on  transpiration  cooling  requirement  is 
conservatively  high  as  the  flow  is  assumed  to  radiate  as  a 
black  body  and  the  effects  of  thermal  radiation  blockage  by  the 
seeded  layer  in  the  core  is  neglected. 

6.8  Secondary  Nitrogen  Flow  ' 

The  heat  flow  through  the  wall  must  be  absorbed  by 
the  flow  in  the  coolant  tubes.   The  flow  properties  in  the 
tubes  are  set  by  matching  the  heat  transfer  at  the  throat 
as  this  is  the  point  of  maximum  heating  load. 

The  flow  in  the  tubes  can  not  be  set  point  by  point 
to  enable  maximum  heat  transfer  to  be  possible  along  the 
entire  length  of  the  nozzle  and  heat  transfer  rates  below 
Q**ool    are  possible  along  the  nozzle.    This  problem  is 
recognized  but  not  investigated  further  in  this  analysis 
as  it  is  felt  that  the  transpiration  cooling  flow  can  be 
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tailored   to  fit  a  wall  with  variations  in  heat  transfer 
capabilities . 

The  heat  flow  through  the  walls  is  absorbed  by  the 
coolant  flow, 

Q™„  -  V,c  (Tw-Tc)  (6.21) 

Here  radiation  and  nuclear  heating  of  the  coolant  are  ignored. 
The  coolant  is  flowing  in  a  nozzle  shaped  duct  but  here 
fluid  is  bled  out  of  the  tube.   In  the  absence  of  data  on 
sucking  of  a  turbulent  boundary  layer,  the  data  of  Reference  11 
is  extrapolated  to  include  the  sucking  case.   This   is  also 
illustrated  in  Figure  10,  and  is  analogous  to  the  treatment 
of  a  sucked  laminar  boundary  layer  given  in  Reference  25. 

The  method  of  Bartz  is  used  as  before  to  determine  the 
heat  transfer  coefficient.   The  coolant  tube  nozzle  throat 
radius  is  set  at  .5  cm. 

The  chamber  temperature  and  pressure  of  the  secondary 
nitrogen  flow  is  set  by  matching  the  heat  transfer  through 
the  wall  at  the  throat  to  the  heat  absorbed  by  the  coolant 
of  that  point,  taking  into  account  the  effect  of  the  transla- 
tion cooling  flow  sucking  the  boundary  layer  in  the  coolant 
tubes. 

For  the  conditions  of  this  study,  the  secondary  nitrogen 
chamber  pressure  is  1010  atmospheres  and  the  temperature  is 

530  K.   With  these  chamber  conditions  the  mass  flow  of 

5 
secondary  nitrogen  is  2.17  x  10   grams  per  second  and  the 

ratio  of  secondary  to  primary  nitrogen  flow  is   0.67  to  1. 
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The  secondary  flow  of  nitrogen  that  is  used  to  provide 
transpiration  and  backside  cooling  of  the  nozzle  is  approx- 
imately the  same  as  the  primary  flow  through  the  nozzle.   This 
doubles  the  amount  of  pump  work  required,  but  as  the  nozzle 
would  fail  without  the  transpiration  protection,  this  additional 
pump  work  is  an  acceptable  penalty  that  must  be  incorporated 
into  the  design  of  the  system. 
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CHAPTER  VII 
AUXILIARY  EQUIPMENT 

7.1  General 

Operating  this  facility  as  a  closed  cycle  removes 
some  of  the  technological  problems  involved  in  the  design 
of  a  gaseous  core  rocket.   The  problem  of  containing  the 
gaseous  fuel  in  the  reactor  core  is  avoided  here  by  operating 
as  a  closed  cycle.   This  introduces  the  problem  of  additional 
equipment  in  the  cycle.   Among  the  additional  equipment  needed 
are  high  powered  turbines,  high  capacity  pumps,  large  heat 
exchangers  and  fuel  separation  and  processing  equipment. 

7.2  Water  Diffusor 

The  flow  at  the  exit  of  the  test  section  is  mixed 
with  the  secondary  nitrogen  flow  which  was  utilized  to  provide 
transpiration   and  backside  cooling  of  the  nozzle.  '  At  this 

point  the  flow  stagnation  temperature,  determined  from  the 

o 

enthalpy  of  the  mixture,  is  8000  K  and  the  stagnation  pressure, 

assuming  isentropic  expansion  through  the  nozzle,  is  1000 
atmospheres .Conventional  heat  exchangers  can  not  operate 
at  these  temperatures  and  pressures,  neither  can  conventional 
turbo-machinery.   To  lower  the  stream's  temperature  and  pressure, 
a  water  spray  diffuser  is  employed.   The  water  is  sprayed 
directly  into  the  hot  nitrogen  stream  and  is  heated  to  super  - 
heated  steam  by  the  nitrogen,  which  loses  energy  in  the  process. 
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The  cross  sectional  area  of  the  diffuser  can  be  designed 
to  obtain  the  desired  outlet  conditions. 

Conservation  of  energy  determines  the  required  water 
flow.   Conservation  of  energy  requires  that 

(7.1) 

where   yy\   is  the  mass  flow  rate 

H   is  the  total  enthalpy 
I    is. the  mole  fraction  of  water 
l~l   is  the  mole  fraction  of  nitrogen 
and  subscript  w   refers  to  water 
i    inlet  conditions 
0         stagnation  conditions 
Z         outlet  conditions. 

Equation  7.1  was  utilized  to  obtain  the  required  mass 
flow  of  water  necessary  to  reduce  the  temperature  and  pressure 
of  the  nitrogen  stream  to  the  desired  values. 

The  size  of  the  water  diffuser  is  set  by  the  perfect 

gas  law,  written  in  the  form 

» 

R   -  _A.  m»  J%     I     Ta  V,  J(,  (7.2) 

where  all  values  are  static  values  and 

c/t   is  the  mixture  average  molecular  weight 
$■       is  the  mixture  average  specific  heat     ratio 

With  the  properties  of  nitrogen  given  in  Reference 
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14     and  those  of  water  from  Reference  26  ,  Equations  7.1 
and  7.2   can  be  used  to  find  the  necessary  water  flow  and 
diffusor  exit  area,  once  the  exit  Mach  number,  temperature 
and  pressure  are  specified. 

The  diffusor  outlet  total  temperature  is  set  at  1000  K, 
as  an  upper  operating  temperature  for  turbo-machinery.   The 
outlet  total  pressure  is  set  a   100  atmospheres  and  the 
Mach  number  at  0.01 

With  these  conditions  set,  the  required  mass  flow  of 

5 
water  is  9  x  10   grams  per  second,  and  the  output  area 

of  the  diffusor  is  8.25  x  10   cm. 

7.3   Fuel  Separator 

The  cooling  process  in  the  water  diffuser  cools  uranium 
fuel  in  the  flow  well  below  its  melting  point.   The  super- 
cooled uranium  rapidly  condenses  and  forms  small  droplets 
and  particles  which  are  suspended  in  the  flow.8  The  thermal 
energy  of  the  flow  is  used  to  provide  a  centrifugal  force 

field  to  separate  the  particles  from  the  gas.  The  drag  on 

25 

these  particles  is  given  for  Stokes  flow  as 

E>  .  4^V  (7.3) 

where  ^a  is  the  viscosity  of  the  steam,  nitrogen  mixture. 
The  equation  of  radial   motion  of  a  uranium  particle  in 
a  centrifugal  force  field  is,  with  the  multiplier  of  velocity 
in  Equation  7.3defined  as  Co, 

r  +C0l    -  *>*r  -  o  (7.4) 

This  has  a  solution  of  the  form 
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r»=5==.<e       2-  *  _  o  2        ~M  (7.5) 


-7\e      *•         S.e         »         'J 


For  small  particles,  the  drag  coefficient  Co  is 
large  and  the  second  exponential  term  of  Equation  7.5 
decays  rapidly  regardless  of  the  value  of  Oi   and  the 
first  term  determines  the  radial   position  of  the 
particle.   For  small  values  of  CO  /the  particle  is  strongly 
restrained  by  the  drag  forces  of  the  fluid  and  a  centrifugal 
separator  is  not  very  effective.  For  a  centifugal  separator 
to  efficiently  separate  the  small  particles,  the  rotational 
velocity  must  be  at  least  of  the  order  of  the  drag  cofficient. 
Under  the  conditions  of  the  flow  just  after  the  water  diffusor, 
the  angular  velocity  must  be  of  the  order  of  100  radians 
per  second  for  micron  sized  particles.   This  imposes  large 
velocity  requirements  if  the  flow  is  swirled  in  one  section 
at   the  outlet  of  the  water  diffusor,  so  the  flow  is  divided 
into  several  centrifugal  separators. 

7.4  Turbines  and  Pumps 

The  ideal  power  available  in  the  flow  as  it  comes  out 
of  the  uranium  separators  is  given  by 

Sf^  YYU  HH  +™~  Hw  (7.6) 

It  is  convenient  however  to  expand  the  flow  only  until 
the  water  vapor  in  the  flow  reaches  its  saturation  point. 
At  this  point,  the  water  can  be  condensed  and  removed  from 
the  flow  stream  and  the  nitrogen  further  cooled  to  reduce 
the  pump  work  required  to  increase  its  pressure  to  core 
entry  conditions.   Under  these  conditions  Equation  7.6 
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becomes 


-*-*%.,  ]  H*  (tocxTk,  |»OOaVml    -  Hw  (  l,AT/  P»AT)   |       (7.7) 


For  isentropic  expansion  of  the  water  vapor,  the  ideal 
power  available  is  2,520  megawatts. 

After  condensation,  the  water  is  pumped  as  a  liquid 

to  120  atmospheres  and  injected  into  the  water  diffusor. 

This  requires  an  ideal  pump  work, 

o 
VPw  -  A£.w  A^  (7.8) 

Ideally,  10.9  megawatts  are  required  for  the  water  pump 
work. 

o 

The  nitrogen  is  cooled  in  a  heat  exchanger  to  293  K 

and  then  pumped  back  to  1000  atmospheres.   This  must  be 
done  in  at  least  two  pumps  with  interpump  cooling  to  maintain 
the  temperature  at  or  below  a  structurally  safe  limit.   The 
ideal  nitrogen  pump  work  is  given  for  each  pump  by 

VHU  =  ■«-lHA.*>H.XM   (7.9) 

where  the  temperatures  and  pressures  are  related  by  the 
isentropic  compression  relation 


T,.    Vt>;  ) 


*  (7.10) 
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The  total  pump  work  for  all  nitrogen  pumps  is  868 
megawatts . 

The  total  heat  rejected  in  the  cycle  is  the  sum 
of  the  heat  of  condensation  of  the  water  vapor,  the 
cooling  of  the  nitrogen  prior  to  being  pumped,  and  inter- 
pump  cooling  after  pumping  so  that  the  nitrogen  can  be 

used  to  cool  the  reactor  moderator.   The  total  heat  rejected 

3 

in  the  cycle  is  7.1  x  10   megawatts. 

If  water  is  used  as  a  heat  sink  for  this  rejected  heat, 
the  mass  flow  required  is 


T« 


i   c^izr 


.        .  o 

For  a  temperature  rise  of  this  coolant  water  of  100  K, 

the  mass  flow  required  is  1.1  x  10   grams  per  second. 


7.5  Minimum  Power  Level 

The  minimum  operating  power  level  of  this  test  facility 
is  set  by  the  criticality  requirements  of  the  reactor,  Chapter 
III,  and  heat  transfer  considerations  in  the  reactor  and  nozzle, 
Chapter  V  and  VI.   This  high  operating  power  requires  high 
mass  flow  rates  and  large  capacity  pumps  and  turbines,  as 
indicated  in  Section  7.4. 

There  is  no  lack  of  power  available  to  drive  this  machinery 
and  small  units  can  operate  in  parallel  to  provide  the  desired 
capacity.   As  the  system  is  proposed,  there  is  a  requirement 
that  the  combined  efficiency  of  the  turbines  and  pumps  be 
about  35  percent,  but  if  all  of  the  energy  in  the  cycle  is 
used,  efficiencies  as  low  as  12  percent  are  acceptable. 
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Equations  5.1  and  6.1  can  be  combined  to  express  the 
operating  power  minimized 'with  respect  to  the  area  of  the 
throat, 

f>  =  (€<rT04  +.o8  IP/A.)  Ac        (7.12) 

The  area  of  the  core  is  minimized  for  any  pressure 
and  temperature  by  setting  the  ratio  of  nitrogen  mass 
flow  to  uranium  mass  flow  at  15.4,  as  was  shown  in  Section  3.2. 

Equation  3.11  shows  that,  with  the  mass  flow  ratio 
fixed,  the  critical  radius,  R,  varies  inversely  with  the 
number  density  of  uranium,  N„   .   For  constant  operating 
pressure,  the  number  density  of  uranium  varies  inversely 
with  the  operating  temperature,  1^   .   The  critical  radius 
is  then  proportional  to  the  operating  temperature  and  the 
core  area  is  proportional  to  the  square  of  the  operating 
temperature. 

With  fixed  mass  flow  ratio  and  operating  pressure,  the 
minimum  operating  power  is,  to  the  first  order, 

P~„  -  C  Ta'  (7.13) 

where  C  is  a  proportionality  constant. 

Equation  7.13  is  valid  only  in  the  temperature  range 
where  the  heat  loads  to  the  moderator  are  too  large  to  be 
removed  by  conventional  cooling  and  the  energy  absorbed  in 
protecting  the  moderator  walls  must  pass  through  the  nozzle. 

A  plot  of  the  estimated  minimum  operating  power  as  a 


68 

function  of  temperature  is  given  in  Figure  11.   At  operating 

o 

temperatures  below  7000  K,  this  facility  is  in  the  power 

category  of  large  arc  tunnels  and  MHD  devices,  and  the  i 

relative  merits  of  each  type  of  facility  would  have  to  be 
considered  before  deciding  which  system  is  most  desirable. 
At  the  higher  temperatures  at  which  this  facility  is  con- 
cept ionally  most  useful,  the  minimum  operating  power  rises 
rapidly,  as  well  the  size  of  the  physical  plant.   It  is  also 
noted  that  reducing  the  operating  temperature  reduces  the 
re-entry  simulation  capability  of  the  test  facility,  as 
illustrated  in  Figure  2. 


The  economics  of  developing  such  a  facility  must  be 
examined  to  determine  if  the  need  for  such  a  testing  capability 
justifies  the  expense  of  so  large  an  undertaking. 
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CHAPTER  VIII 
CONCLUSIONS 

Within  the  scope  of  the  material  examined  in  this 
report  it  appears  to  be  within  the  capability  of  projected 
technology  to  design  and  build  a  high  enthalpy  test  facility 
powered  by  a  gaseous  core  reactor. 

This  test  facility  can  not  exactly  duplicate  all  of 
the  conditions  imposed  upon  a  vehicle  during  an  atmospheric 
re-entry,  but  does  provide  a  means  of  testing  re-entry 
materials  and  shapes  under  the  high  temperature  and  shear 
conditions  which  are  encountered  during  a  re-entry.   Existing 
facilities  cannot  attain  the  high  temperature,  high  pressure, 
long  duration  flows  necessary  to  test  the  shear  mode  of 
failure.   The  gaseous  core  reactor  provides  the  power  required 
for  these  flows. 

Nuclear  contamination  of  several  proposed  re-entry 
materials  is  slight  enough  so  that  the  model  may  be  examined 
immediately  after  a  test  run.   This  may  not  be  true  in 
general  and  each  proposed  test  material  must  be  examined 
by  the  methods  of  Chapter  IV  to  determine  the  degree  of 
contamination . 

The  high  neutron  absorption  cross  section  of  nitrogen 
as  a  working  fluid  causes  the  reactor  to  be  large  and  the 
mass  flow  ratio   of  uranium  to  nitrogen  to  be  relatively 
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high-   The  reactor  size  and  operating  power  is  minimized 
by  setting  the  mass  ratio  of  nitrogen  to  uranium  in  the 
core  at  15.4  to  1.   This  adversely  affects  the  quality 
of  the  mixture  for  simulating  atmospheric  conditions, 
but  is  slight  enough  to  be  considered  acceptable. 

The  minimum  operating  power  of  the  test  facility  is 
set  by  the  radiant  power  of  the  core  which  cannot  be 
absorbed  by  the  moderator  and  must  pass  through  the  nozzle. 
In  the  neighborhood  of  10,000   K,  for  a  constant  mass  ratio, 
the  minimum  operating  power  varies  as  the  temperature  in 
the  core  to  the  sixth  power  and  as  the  inverse  of  the  core 
pressure . 

Operation  of  the  facility  at  the  minimum  power  level 
requires  high  mass  flow  rates  and  associated  high  capacity 
pumps  and  turbines.   These  pumps  must  be  capable  of  high 
pressure  and  high  pressure  ratio  operation.   Efficiency 
of  the  combined  turbine  and  pump  can  be  as  low  as  twelve 
percent  for  self  sustaining  operation  of  the  system. 

The  turbo-machinery  required  set  a  practical  upper 
limit  on  the  operating  power  of  the  system.  Structural 
considerations  impose  an  upper  pressure  limit  on  the  system. 

This  analysis  has  examined  only  the  minimum  powered 
point  for  the  operation  at  the  maximum  pressure  and  temperature 
and  has  found  that  such  a  facility  appears  technically  feasible 
Operation  of  a  facility  at  lower  power  levels  is  possible 
with  a  reduction  in  chamber  temperature  and  associated  loss 
of  re-entry  simulation  capability. 
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It  is  recommended  that  the  economics' involved  in  the 
construction,  maintenance  and  operation  of  this  system  be 
examined  to  determine  if  such  a  large  facility  would  be 
economically  practical. 
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TABLE  I 


FISSION  PRODUCTS 


Radionuclide   Half  Life   Yield  Fraction   Decay  Constant 


A 

£*»(■«  ) 

* 

AACseO 

Br 

87 

55.6 

io"3 

1.2  x  IO"2 

Br 

89 

4.51 

5  x  10"3 

1.54  x  IO"1 

Ag 

110 

24.2 

io"4 

2.9  x  IO"2 

In 

114 

72 

io"4 

9.6  x  IO"3 

In 

116 

13 

io"4 

5.4  x  10~2 

Sb 

135 

1.52 

4  x  IO"2 

4.56  x  IO"1 

I 

137 

22 

6  x  10~2 

3.1  x  IO"2 
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